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Abstract. Using recently evaluated contributions (including a novel one calculated here), we present 
updated values for the pole mass and MS mass of the b quark: mi, = 5022 ± 58 MeV, for the pole 



o 



43 



mass, and mj,(mj) = 4286±36 MeV for the MS one. These values are accurate including, respectively, 
0(a 5 s loga s ) and 0(al) corrections and, in both cases, leading orders in the ratio m\l'm%. 



^C) ' One of the sources of information for the quark masses is quarkonium spectroscopy. By evaluating the 

OA • bb potential including relativistic and radiative corrections, as well as leading nonperturbative effects, I 1,2 ' 3 ' 

and using this in a perturbative expansion, it has been possible to find values of the pole quark masses with 
increasing accuracy;' 2,3 ' 4 ' in this note we will go up to fourth and leading fifth order, in the approximation of 
' neglecting "light" (u, d, s) quark masses and to leading order (actually, 0{a s m 2 c /m 2 )) in the c quark mass. 

£T) . The connection with the MS mass has been known for some time to one and two loops' 5 ! : very recently, a 

three loop evaluation has been completed. Coupling this with the pole mass evaluations, we now have an 
order a 3 result for the MS mass. We review here briefly this. 

o 
o 
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1. nib — m(,(mb) connection 
Write, for a heavy quark, 



p m{m) =m/{l + 8 1 + 6 2 + 6 a -\ }; (la) 

O i' m here denotes the pole mass, and m is the MS one. One has 

<D ' 

— * = C«!M, * 2 = c 2 (^) 2 , <^c 3 (^) 3 . (lb) 

^ ' Here a s is to be calculated to three loops: 

ft ' 

. 4tt f AlogL /? 2 log 2 L-/? 2 logL + /3 2 /3 -/3 2 
aAfi) = M\ 1 - I* 

with 



(3 L { faL ffiL 2 

u 2 

T — l^rr ^ . a _ 11 2„ a _ i n o 38„ a _ 2857 5033 „ , 325„2 

L - l0 %^2' P0-H-3"-/, ft -102 --3- n f , P2-— i8~ n / + "54 n /- 



A 2 

The coefficient C2 has been evaluated by Gray et &lS 5 \ and reads 



c 2 = -K + 2C F , (lc) 



K =Ko + £ A O > K o = ^ 2 lo g 2 + ^ 2 - I C (3) + W ~ for 3 + &) (»/ + 1) 

i=l 

-16.11- 1.04 n f ; A(p) = § [^p - \p 2 + • ■ ■] . 



(Id) 



m,i are the (pole) masses of the quarks strictly lighter than m, and n/ is the number of these. For the b 
quark case, n/ = 4 and only the c quark mass has to be considered; we will take m c = 1.8 GeV (see Table 1 
below) for the calculations. 



1 



8i(b) 


=0.090, 


81(c) 


=0.137, 


62(b) 


=0.045, 


82(c) 


=0.108, 


83(b) 


=0.029; 


83(c) 


=0.125. 



The coefficient C3 was recently calculated by Melnikov and van Ritbcrgen,! 6 ! where the exact expression 
may be found. Neglecting now the rrii, 

c 3 ~ 190.389 - 26.655171/ + 0.652694n). (le) 

For the b, c quarks, with a s as given below, 

.090, 5i(c) =0.137, 

(2) 

From these values we conclude that, for the c quark, the series has started to diverge at second order, 
and it certainly diverges at order a 3 ,. For the b quark the series is at the edge of convergence for the a 3 
contribution. 

Take now as input parameters 

A(n f = 4, three loops) = 0.283 ± 0.035 GcV [ a s (M§) ~ 0.117 ± 0.024 ] 

(ref. 7) and for the gluon condensate, very poorly known, the value (a s G 2 ) = 0.06 ± 0.02 GeV 4 . From the 
mass of the Y particle we have a very precise determination for the pole mass of the b quark. This determina- 
tion is correct to order and including leading 0(m 2 /m 2 ) and leading nonperturbative corrections as well 
as the a 5 s corrections proportional to loga s ; the details of it will be given below. With the renormalization 
point fj, — rahCpds we have, 

m b =5022 ± 43 (A) =F 5 ({a s G 2 ))+ll (vary ^ by 25%) ± 38 (other th. unccrt.) , 

(3a) 

=5022 ± 58 MeV . 

Here we append (A) to the error induced by that of A, and likewise ((a s G 2 )) tags the error due to that of 
the condensate. The error labeled (other th. uncert.) includes also the error evaluated in ref. 8; the rest is 
as in ref. 3. 

Using the three loop relation (1) of the pole mass to the MS mass we then find 

m b (m b ) = 4284 ± 7 (A) =p 5 ({a s G 2 }) ± 35 (other th. unccrt.) = 4284 ± 36 MeV . (3b) 

The slight dependence of m on A when evaluated in this way was already noted in ref. 2. 

There is another way of obtaining m, which is to express directly the mass of the T in terms of it, using 
Eq. (1) and the order a 3 formula for the T mass in terms of the pole mass (see e.g. ref. 2). One finds, for 
rif =4, and neglecting m 2 /m 2 , 

M(T) = 2m(m) <jl + C F ^^ + 7.559 f-MV + [66.769+ 18.277 (log C F + log a.(m))\ ( 

J 7T \ 7T J \ 7T 

(4a) 

(One could add the leading nonperturbative contributions to (4a) a la Leutwyler-Voloshin in the standard 
way; see e.g. refs. 2, 3, 9). This method has been at times advertised as improving the convergence, allegedly 
because the MS mass does not suffer from nearby renormalon singularities. But a close look to (4a) does not 
seem to bear this out. To an acceptable 0(af) error we can replace \og(a s (m)) by log(a s (M(T/2)) above. 
With A as before (4a) then becomes 

a s (fh) ( a s (fh,y 2 



M(T) = 2m 6 (m fe )|l + C F ^^ + 7.559 (^f^-J + 43.502 (^) j. (4b) 

This does not look particularly convergent, and is certainly not an improvement over the expression using 
the pole mass, where one has for the choice' 3 ! [i = CFin b a s , and still neglecting the masses of quarks lighter 
than the b, 

M(T) = 2m b j 1 - 2.193 (^^) ' - 24.725 (^^) ' - 458.28 (^^) ^ + 897.93 [loga s ] (^)" j . (5) 
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To order three, (5) is actually better 1 than (4b). What is more, logarithmic terms appear in (4) at order a^, 
while for the pole mass expression they first show up at ajj. Finally, the direct formula for M(T) in terms of 
the MS mass presents the extra difficulty that the nonperturbative contribution becomes larger than than 
what one has for the expression in terms of the pole mass (~ 80 against ~ 9 MeV) , because of the definition of 
the renormalization point. With the purely perturbative expression (4) plus leading nonperturbative (gluon 
condensate) correction one finds the value fhb(fhb) = 4167 MeV, rather low. 

2. Improved determination of m b 

Eq. (5) was deduced neglecting the masses of all quarks lighter than the b. The influence of the nonzero 
mass of the c quark, the only worth considering, will be evaluated now. To leading order it only contributes to 
the bb potential through a c-quark loop in the gluon exchange diagram (diagram fi in ref. 2) . The momentum 
space potential generated by a nonzero mass quark through this mechanism is then, in the nonrelativistic 
limit, 

v c _ - -ggffii /' - ,) ,o 6 + * ~ ^ . (6) 

k Jo A* 
We expand in powers of ml/Is 2 . The zeroth term is already included in (5). The first order correction is 

~ _ 8C F T F a 2 s m 2 c 

Oc mass v — .^4 ' V / 



In x- space, 



6 emm V = C * T ** m ° r. (8) 



This induces the shift in the mass of the T of 



^cmassM(r) = § m bl 



so Eq. (5) is modified to 



WTP. =2m b { 1 - 2.193 (^M.) - 24.725 ( ^iMV - 458.28 ( " J '" 



+897.93 [loga s ] (^) 



7T 

5 3T F a s m\ 



(8) 



2tt ml 



This produces the value quoted in (3a). Note that the (new) correction of order m 2 /m 2 is responsible for a 
shift in TOfc of 

^cmass^b = -35 MeV, 

substantially larger than the ct^ log a s correction evaluated by Brambilla et al. M which, for the renormaliza- 
tion point /i = mbCpOis, gives 

S[alio gas ]m b = \m b [C F + %CA\C F a 5 s (loga s )/ir ~ -8 MeV. 

We collect in the table the determinations of the b quark mass based on spectroscopy, to increasing 
accuracy. The stability of the numerical values of the pole mass is remarkable: the pole masses all lie within 
each other error bars. The MS ones show more spread. 



1 The convergence of Eq. (5) is still improved if one solves exactly the purely coulombic part of the static potential, 
as was done in refs. 2, 3, where we send for details. For example, the O(af) term becomes — 232.12(a s /7r) 4 . This 
is the method we used to get the values of m b here. 
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Reference 


TOb(pole) 


m b {ml) 


m c (pole) 


m c (rr? c ) 


TY 
PY 
Here 


4971 ± 72 
5065 ± 60 
5022 ± 58 


4401^35 
4455±g 
4286 ± 36 


1585±20(*) 
1866t£jj 


1321±30(*) 
1542±ig 



Table 1. b and c quark masses. (*) Systematic errors not included. 



TY: Titard and Yndurain [2] . 0(a 3 s ) plus 0(a 3 s )v, 0(v 2 ) form; 0(a 2 B ) form. Rescaled for A{n f = 4) = 283 MeV. 
PY: Pineda and Yndurain 131 . Full O(aj) for m; 0(a 2 s ) for fa. Rescaled for A(n f = 4) = 283 MeV. 
Here: This calculation. O(af), 0(a s mc/mj) and 0(Qsloga s ) for m; 0(a?,) and O(o^iric/m[) for fa. Values 
not given for the c quark, as the higher order terms are as large as the leading ones. 



We finally remark that the values of m& quoted e.g. in the Table 1 were not obtained solving Eq. (8), but 
solving exactly the coulombic part of the interaction, and perturbing the result (see refs. 2, 3 for details). 
We also note that, in the determinations of m&, the new pieces, 0(a s m^/mf) and 0(a^ loga s ), have been 
evaluated to first order; in particular, we have included the corresponding shifts in the central values, not in 
the errors. If we included these, the errors would decrease by some 7%. 
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